Quantum Monte Carlo study on speckle variation due to photorelaxation of 
ferroelectric clusters in paraelectric barium titanate 
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Time-dependent speckle pattern of paraelectric barium titanate observed in a soft x-ray laser 
pump-probe measurement is theoretically investigated as a correlated optical response to the pump 
and probe pulses. The scattering probability is calculated based on a model with coupled soft x-ray 
photon and ferroelectric phonon mode. It is found that the speckle variation is related with the 
relaxation dynamics of ferroelectric clusters created by the pump pulse. Additionally, critical slowing 
down of cluster relaxation arises on decreasing temperature towards the paraelectric-ferroelectric 
transition temperature. Relation between critical slowing down, local dipole fluctuation and crystal 
structure are revealed by quantum Monte Carlo simulation. 

PACS numbers: 78.47.-p, 61.20.Lc, 63.70.+h, 77.80.-c 
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I. INTRODUCTION 

Speckle is the random granular pattern produced when 
a coherent light is scattered off a rough surface. It car- 
ries information of the specimen surface, for the intensity 
and contrast of the speckle image vary with the rough- 
ness of surface being illuminated. 1 Numerous approaches 
have been devised to identify surface profiles by either 
the speckle contrast or the speckle correlation method. 2 
Recent application of pulsed soft x-ray laser has improved 
the temporal and spatial resolution to a scale of picosec- 
ond and nanometer. By this means, dynamics of surface 
polarization clusters of barium titanate (BaTiOa) across 
the Curie temperature (T c ) has been observed; 3 ^ which 
paves a new way to study the paraelectric-ferroelectric 
phase transition. 

As a prototype of the ferroelectric perovskite com- 
pounds, BaTi03 undergoes a transition from paraelectric 
cubic to ferroelectric tetragonal phase at T c =395 K. Be- 
low T c , two kinds of ferroelectric domain are developed 
with mutually perpendicular polarization. Structural 
phase transition and domain induced surface corrugation 
have been observed via atomic force microscopy, 5 scan- 
ning probe microscopy^ neutron scattering^ and polar- 
izing optical microscopy- In addition to the extensive ap- 
plication of BaTiC>3 in technology due to its high dielec- 
tric constant and switchable spontaneous polarization, 9 
there is also an enduring interest in understanding the 
mechanism of paraelectric-ferroelectric phase transition. 
It is generally considered that the transition is a classic 
displacive soft-mode type driven by the anharmonic lat- 
tice dynamics! 10 ' 11 However, recent studies have also sug- 
gested an order-disorder instability which coexists with 
the displacive transition ! 12 ' 13 Therefore, direct observa- 
tion on creation and evolution of ferroelectric cluster 
around T c is of crucial importance for clarifying the na- 



ture of phase transition. Since the above-mentioned con- 
ventional time-average-based measurements cannot be 
adapted to the detection on ultrafast transient status 
of dipole clusters, diffraction speckle pattern of BaTiC>3 
crystal measured by the picosecond soft x-ray laser has 
turned out to be an efficient way for this purpose. 

Very recently, Namikawa et al^ study the polarization 
clusters in BaTiC>3 at above T c by the plasma-based x-ray 
laser speckle measurement in combination with the tech- 
nique of pump probe spectroscopy. In this experiment, 
two consecutive soft x-ray laser pulses with wavelength of 
160 A and an adjustable time difference are generated co- 
herently by the Michelson type beam splitter. After the 
photo excitation by the pump pulse, ferroelectric clusters 
of nano scale are created in the paraelectric BaTiC>3 and 
tends to be smeared out gradually on the way back to the 
equilibrium paraelectric state. This relaxation of cluster 
thus can be reflected in the variation of speckle intensity 
of the probe pulse as a function of its delay time from the 
first pulse. It has been found that the intensity of speckle 
pattern decays as the delay time increases. Moreover, 
the decay rate also decreases upon approaching T c , in- 
dicating a critical slowing down of the cluster relaxation 
time. Hence, by measuring the correlation between two 
soft x-ray laser pulses, the real time relaxation dynam- 
ics of polarization clusters in BaTiC>3 is clearly repre- 
sented. In comparison with other time-resolved spectro- 
scopic study on BaTi03, for example the photon correla- 
tion spectroscopy with visible laser beanifi 5 - Namikawa's 
experiment employs pulsed soft x-ray laser as the light 
source. For this sake, the size of photo-created ferroelec- 
tric cluster is reduced down to a few nanometers, and the 
cluster relaxation time is at a scale of picosecond. This 
measurement, thus, offers a new insight into the ultrafast 
quantum dynamics of domain structure. 

In this work, we examine the above-mentioned novel 
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behaviors of ferroelectric cluster observed by Namikawa 
from a theoretical point of view, aiming to provide a basis 
for understanding the critical nature of BaTi03. Theo- 
retically, the dynamics of a system can be adequately de- 
scribed by the linear response theory, i.e., to express the 
dynamic quantities in terms of time correlation functions 
of the corresponding dynamic operators. In general, the 
path integral quantum Monte Carlo method is computa- 
tionally feasible to handle the quantum many body prob- 
lems, for it allows the system to be treated without mak- 
ing any approximation. However, simulation on real time 
dynamics with Monte Carlo method is still an open prob- 
lem in computational physics because of the formidable 
numerical cost of path summation which grows expo- 
nentially with the propagation time. The common ap- 
proach to circumvent this problem is to perform imagi- 
nary time path integration followed by analytic continu- 
ation, and to compute the real time dynamic quantities 
using Fourier transformation. In the present study, the 
real time correlation functions and real time dependence 
of speckle pattern are investigated by this scheme. Our 
quantum Monte Carlo simulation demonstrates that the 
relaxation dynamics of photo-created nano clusters plays 
an essential role in determining the delay time depen- 
dence of speckle variation. Furthermore, it is found that 
the critical slowing down of photorelaxation is related 
to the local dipole fluctuation, which arises near T c and 
stablizes the photo-created ferroelectric cluster. 

The remaining of the present paper is organized as fol- 
lows. In Sec. II, the model Hamiltonian and theoretical 
treatment are elaborated. In Sec. Ill, our numerical re- 
sults on speckle correlation and critical slowing down are 
discussed in details. In Sec. IV, a summary with conclu- 
sion is presented finally. 

II. THEORETICAL MODEL AND METHODS 
A. Model Hamiltonian 

The theoretical interpretations for structural phase 
transition and domain wall dynamics have be well estab- 
lished in the framework of Krumhansl-Schrieffer model 
(also known as <fi 4 model) i 16 ' 17 i 18 i 19 In this model, the 
particles are subject to anharmonic on-site potentials and 
harmonic inter-site couplings. The on-site potential is 
represented as a polynomial form of the order parameter 
such as polarization, displacement, or elasticity, which 
displays a substantial change around T c . Since the <j) 
model is only limited to second-order transitions, in the 
present work we invoke a modified Krumhansl-Schrieffer 
model (also called </> 6 model) 20 ' 21 to study the first-order 
ferroelectric phase transition of BaTiOa. In this scenario, 
the Hamiltonian of BaTiOa crystal (= Hf) is written as 
(here we let % = 1), 



u ° = y £(<??- <*Q? + f<??). (2) 

(1,1') 

where, Uq and Ui are the on-site potential and inter-site 
correlation, respectively. Qi is the coordinate operator 
for the electric dipole moment due to a shift of titanium 
ions against oxygen ions, i.e., the T\ u transverse optical 
phonon mode, ojq is the dipole oscillatory frequency, I la- 
bels the site, and (I, I') in Eq. (3) enumerates the nearest 
neighboring pairs. 

In order to describe the optical response of BaTi03 
due to x-ray scattering, we design a theoretical model 
to incorporate the radiation field and a weak interplay 
between radiation and crystal. The total Hamiltonian 
reads, 

H = H p + Hf + H p f, (4) 

where 

Hp = y^OfcQfcfflfc, £lk = c\k\, (5) 
k 

is the Hamiltonian of polarized light field, at (dfc) is the 
creation (annihilation) operator of a photon with a wave 
number k and an energy Qk- c is the light velocity in 
vacuum. In Namikawa's experiment, the wave length of 
x-ray is 160 A, thus the photon energy is about 80 eV. 
Denoting the odd parity of T\ u mode, the photon-phonon 
scattering is of a bi-lincar Raman type, 

V \ - + 

H pf = n 1^ a k+z a k-§Q q >-lQ- q '-^ ( 6 ) 

q,q',k 

where V is the photon-phonon coupling strength, Q q (= 
^ e~ %ql Qi) the Fourier component of Qi with a 
wave number q. Without losing generalitivity, here we 
use a simple cubic lattice, and the total number of lattice 
site is N. 



B. Optical response to pump and probe photons 

Since there are two photons involved in the scattering, 
the photon-phonon scattering probability can be written 

as, 

p{t) = E <k> (o)4, ( a h (<)4i ( A + *) 

xa K (A + t)al o (t)a kl (A)al o (0))), (7) 

where 

((•••))=Tr(e-^---)/(e-^), (8) 
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FIG. 1: Pulse sequence in an x-ray laser speckle experiment. 
The pump and probe pulses of ko creates and detects ferro- 
electric clusters in the sample of paraelectric BaTiOa , respec- 
tively, and generate new x-ray fields in the direction ki and 
k[ after a short time interval A. 



means the expectation, f3 (= 1/ksT) is the inverse tem- 
perature, and the time dependent operators are defined 
in the Heisenberg representation, 



0(t) = 



(9) 



Here, t denotes the time difference between two incident 
laser pulses as manifested in Fig. 1, and kg the wave 
number of incoming photon. After a small time interval 
A, the photon is scattered out of the crystal. k\ and k[ 
are the wave numbers of the first and second outgoing 
photons, respectively. 

Treating H p f as a perturbation, we separate Hamilto- 
nian of Eq. (4) as, 



H = Hn 



H 



where 



H — H v + H 



/- 



(10) 



(11) 



is treated as the unperturbed Hamiltonian. By expand- 
ing the time evolution operator in Eq. (9) with respect 
to H p j, 
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(12) 



we find that the lowest order terms which directly depend 
on t arc of fourth order, 

/•A pA pA pA 

P(t) dh dt 2 / dt\ \ dt' 2 V 
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(13) 



where the operators with carets are defined in the inter- 
action representation, 



0(t) 
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(14) 



Fig. 2 represents a diagram analysis for this phonon- 
coupled scattering process, where photons (phonons) are 
depicted by the wavy (dashed) lines, and the upper 
(lower) horizontal time lines are corresponding to the bra 
(ket) vectors^ Diagram (a) illustrates the changes of 
wave number and energy of photons due to the emitted 
or absorbed phonons. This is noting but the Stokes and 
anti-Stokes Raman scattering. Whereas, diagrams (b)- 
(e) are corresponding to the exchange, side band, rapid 
damping and rapid exchange effects, respectively. 

Obviously, diagram (c) brings no time dependence, 
while diagrams (d) and (e) only contributes a rapid re- 
duction to the time correlation of two laser pulses because 
of the duality in phonon interchange. In this sense, the 
time dependence is primarily determined by the diagrams 
(a) and (b). Thus, the scattering probability turns out 
to be, 
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X {{Qq> {A)Q-q+q' (* + 4) 

xQ_ q , (t + t' 2 )Q q , (t + t 2 )Q q - q > (t + h) 

xQ- q+q <(ti)Q- q ,(h))), (15) 

where the photons and phonons are decoupled, and it 
becomes evident that the origin of the i-dependence is 
nothing but the phonon (dipole) correlation. 

Since the photonic part in Eq.(15) is actually time- 
independent, and in the case of forward x-ray scattering 
we have |fco|~|fci|~|^|, the normalized probability can 
be simplified as, 



p(t) _E qq '\((Q 2 q ))G q+q >(t)\' 
P(P) E q , q ,\((Q 2 q Q 2 q+q ,))\ 2 



where 



GM) 



-i2((TQ q (t)QJ0))), 



(16) 



(17) 



is the real time Green's function of phonon, and T the 
time ordering operator. In deriving Eq. (16), we have 
also made use of the fact that the light propagation time 
in the crystal is rather short. The Fourier component of 
Green's function, 



dtG q (t)e~ 



(18) 



is related to the phonon spectral function [= A 9 (o;)] 
through,— 



G fl (w) 



duj' A q (w') 
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(19) 



The phonon spectral function describes the response of 
lattice to the external perturbation, yielding profound 
information about dynamic properties of the crystal un- 
der investigation. Once we get the spectral function, 
the scattering probability and correlation function can 
be readily derived. 



C. Dynamics of crystal 

A mathematically tractable approach to spectral func- 
tion A q (uj) is to introduce an imaginary time phonon 
Green's function, for it can be evaluated more easily than 
its real time counterpart. In the real space, the imaginary 
time Green's function is defined as, 



G u »(t) = -2((TQ 1 {t)Q v ,(P))) ) 



(20) 



where r (= it) is the argument for imaginary time (in 
this paper, we follow a convention of using Roman t for 
real time and Greek r for imaginary time). The imagi- 
nary time dependence of an operator in the interaction 
representation is given by 



0{t) = e THo Oe 



(21) 



Under the weak coupling approximation, and by using 
the Suzuki- Trotter identity, the Green's function can 
be rewritten into a path integral form (here we assume 

G u »(t)= [T>xe-n*'W-*A[-2x,(T)z u ,(0)], (22) 



(23) 



where xi is the eigenvalue of Qi, 

Qi\xi) = x t \xi), 
&f(x) is the path-dependent phonon free energy 

e -/8*/(x) = e - Jo (Ml, [;c(t)] 

with 



(24) 
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(25) 



(l,V) 



and $/ is the total phonon free energy, 



(26) 



In the path integral notation, the internal energy of crys- 
tal Ef (= ((Hf))) is represented as 



E f = 



p xe -^[*/(a!)-*/] 



U 



E 



3 4 
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(27) 



from which the heat capacity can be derived as 



dEf\ 



(28) 



The Green's function of momentum space is given by, 
G q (T) = ±J2 eiqil ~ l " )G »"^ ( 29 ) 



which is connected with the phonon spectral function 
A q {uj) through^ 



G q {r) 



du> cosh — r) w] 
2^ sinh(i/W) 



(30) 



Solving this integral equation is a notoriously ill-posed 
numerical problem because of the highly singular nature 
of the kernel. In order to analytically continue the imag- 
inary time data into real frequency information, special- 
ized methods are developed, such as maximum entropy 
method 2 ^ and least squares fitting method^ In present 
work, we adopt the iterative fitting approach, 24 for it can 
give a rapid and stable convergence of the spectrum with- 
out using any prior knowledge or artificial parameter. 
Since the phonon spectral function does not yield a spe- 
cific sum rule like the case of electron, here we introduce 
an auxiliary spectral function A q (uj) which is defined by 



A g (uj) 



coth(i/3w) 
G q (f3) 



A q (u). 



(31) 



Substituting A q (ui) into Eq. (30), we get 



G q (r) 



duo cosh [(i/? - r) w] 
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cosh (|/3w) 



G,(/3)i». (32) 



It can be easily shown that this auxiliary spectral func- 
tion satisfies a sum rule, 



Z7T 



(33) 



which allow us to solve the integral equation of Eq. (32) 
by the iterative fitting approach. Once A q (ui) is repro- 
duced, the phonon spectral function A q (u>) can be ob- 
tained from Eq. (31). 




FIG. 2: Double-sided Feynman diagrams for scattering process of photon with electric dipole moment (phonon). The photons 
and phonons are denoted by the wavy and dashed lines, respectively. In each diagram, the upper and lower horizontal time 
lines represent the bra and ket vectors, respectively. 



III. NUMERICAL RESULTS AND DISCUSSION 



A. Optical responses 



Based on the path integral formalisms, the imaginary 
time Green's function can be readily calculated via a 
standard quantum Monte Carlo simulation. 24 Our nu- 
merical calculation is performed on a 10x10x10 cubic 
lattice with a periodic boundary condition. The imag- 
inary time is discretized into 10-20 infinitesimal slices. 
As already noticed for the analytic continuation^ if the 
imaginary time Green's function is noisy, the uncertainty 
involved in the inverse transform might be very large, and 
the spectral function cannot be determined uniquely. In 
order to obtain accurate data from quantum Monte Carlo 
simulation, a hybrid algorithm 2 ^ has been implemented 
in our calculation. Besides, we pick out each Monte Carlo 
sample after 100-200 steps to reduce the correlation be- 
tween adjacent configurations. The Monte Carlo data 
are divided into 5-10 sets, from which the 95% confi- 
dence interval is estimated through 10,000 resampled set 
averages by the percentile bootstrap method. We found 
that about 1,000,000 Monte Carlo configurations are suf- 
ficient to get well converged spectral functions and real 
time dynamic quantities. 

In the numerical calculation, the phonon frequency ujq 
is assumed to be 20 meV^ the inter-site coupling con- 
stant g?2 is fixed at a value of 0.032, whereas C4 and c% are 
selected to make the on-site Uq a symmetric triple-well 
potential. As shown in Fig. 3, this triple- well structure 
is featured by five potential extrema located at x a , ±x 




FIG. 3: On-site potential Uo for the modified Krumhansl- 
Schrieffer model in the unit of ujq. The coordinates of poten- 
tial extrema are denoted by x a , %b, and x c . 5 and £ are two 
parameters adopted to characterize this potential. 



and ±x c , where 



x a = 0, 



C4 - \/cl - c 6 



(-6 



C4 + \/c 2 - C 6 



Cf, 



(34) 
(35) 

(36) 



In Figs. 4 and 5, we show the optical responses 
of crystal, where c 4 =2.0132x 10~ 2 and c 6 =3.2595x 10~ 4 
are used. Fig. 4 presents the phonon spectral func- 
tions in the paraelectric phase at different temperatures: 
(a) T=1.001T C , (b) T=1.012T C , (c) T=1.059T C and (d) 
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FIG. 4: Phonon spectral function along the line FXMR of 
Brillouin zone in the paraelectric phase at various tempera- 
tures: (a) l.OOlTc, (b) 1.012T C , (c) 1.059T C , and (d) 1.176T CJ 
where T c =386 K. The inset of panel (a) shows the Brillouin 
zone with high symmetry lines. The inset of panel (b) repre- 
sents the spectrum at F point when T—1.012T C . Error bars 
mark the 95% confidence interval. 



T=1.176T C , where T c =386 K. In each panel, the spectra 
are arranged with wave vectors along the TXMR direc- 
tion of Brillouin zone [see in the inset of panel (a)], and lj 
refers to energy. In the inset of panel (b), the spectrum at 
r point for T=1.012T C is plotted with 95% confidence in- 
terval illustrated by the error bars. Since the spectra are 
symmetric with respect to the origin u>=0, here we only 
show the positive part of them. In Fig. 4, when the tem- 
perature decreases towards T c , as already well-known for 
the displacive type phase transition, the energy of phonon 
peak is gradually softened. In addition, a so-called cen- 
tral peak, corresponding to the low energy excitation of 
ferroelectric cluster, appears at the T point. The collec- 
tive excitation represented by this sharp resonant peak is 
nothing but the photo-created ferroelectric cluster. On 
decreasing temperature, spontaneous polarization is de- 
veloped locally as a dipole fluctuation in the paraelectric 
phase. This fluctuation can stabilize the photo-created 
ferroelectric cluster, leading to a dramatically enhanced 
peak intensity near T c . 

The appearance of sharp peak at T point nearby T c 
signifies a long life-time of the photo-created ferroelectric 
clusters after irradiation. Thus, near T c , they are more 
likely to be probed by subsequent laser pulse, resulting 
in a high intensity of speckle pattern. Keeping this in 
mind, we move on to the results of scattering probability. 
In Fig. 5, we show the variation of normalized probabil- 
ity P(t)/P(0) as a function of t (time interval between 
the pump and probe photons). Temperatures for these 
curves correspond to those in the panels (a)-(d) of Fig. 4, 
respectively. In this figure, P(t)/P(0) declines exponen- 
tially, showing that the speckle correlation decreases with 
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FIG. 5: Normalized speckle scattering probability as a func- 
tion of time for paraelectric BaTiOa , at various temperatures. 
Horizontal dashed line denotes P(t) — P(0)/e. 

t increases as a result of the photorelaxation of ferroelec- 
tric cluster. When t is long enough, the crystal returns to 
the equilibrium paraelectric state. In addition, as shown 
in the figure, the relaxation rate bears a temperature de- 
pendence. On approaching T c , the duration for return 
is prolonged, indicative of a critical slowing down of the 
relaxation. This is because with the decrease of temper- 
ature, the fluctuation of local polarization is enhanced, 
and a long range correlation between dipole moments is 
to be established as well, making the relaxation of photo- 
created clusters slower and slower. 



B. Critical slowing down of photorelaxation 

In order to quantitatively depict the critical slowing 
down, we introduce a relaxation time t r to estimate the 
time scale of relaxation, which is the time for P(t) to 
be reduced by a factor of e from P(0). In Fig. 5, 
P(t)=P(0)/e is plotted by a horizontal dashed line. Cor- 
respondingly, t r is the abscissa of the intersection point 
of relaxation curve and this dashed line. In Fig. 6, the 
relaxation time for various local potential Uq is presented 
at T > T c . Here we adopt two legible parameters, 5 and 
£, to describe the potential wells and barriers for Uq (see 
Fig. 3), which are defined by 

S = [U (x b )-U (x c )]/cj , (37) 
£ = U {x c )/uq. (38) 

Provided 5 and £, C4 and eg can be derived in terms of 
Eqs. (35)-(38). The values of C4 and cq for the calculation 
of Fig. 6 are listed in Table I, where we set £=3.061 and 
change 5 from 4.239 to 4.639. The leftmost point on 
each curve denotes the t r at just above T c , which is a 
temperature determined from the singular point of Cj 
according to Eq. (28). 

As revealed by the NMR experiment^ the 
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FIG. 6: Temperature dependence of relaxation time t r for 
various 8 when T > T c , where £ is fixed at 4.439. Error bars 
show the 95% confidence interval. 



FIG. 7: Temperature dependence of relaxation time t r for 
various £ when T > T c , where 8 is fixed at 3.061. Error bars 
show the 95% confidence interval. 



TABLE I: Parameters adopted for calculation of Fig. 6. 



TABLE II: Parameters adopted for calculation of Fig. 7. 



c 4 


C6 


8 




To (K) 


c 4 


C6 


8 




T c (K) 


2.0696 xl(T 2 


3.4521 xl0~ 4 


4.239 


3.061 


340 


1.9626xl0~ 2 


3.1070xl0~ 4 


4.439 


3.261 


354 


2.0132xl0~ 2 


3.2593xl0~ 4 


4.439 


3.061 


386 


2.0132xl0~ 2 


3.2593xl0~ 4 


4.439 


3.061 


386 


1.9596xl0~ 2 


3.0814xl0~ 4 


4.639 


3.061 


422 


2.0663xl0~ 2 


3.4223xl0~ 4 


4.439 


2.861 


404 



paraelectric-ferroelectric phase transition of BaTiC>3 
has both displacive and order-disorder components in 
its mechanism. Short range dipole fluctuation arises 
in the paraelectric phase near T c as a precursor of the 
order-disorder transition, and condenses into long range 
ferroelectric ordering below T c . Thus, in the present 
study, the relaxation of photo-created cluster is also 
subject to the dynamics of this dipole fluctuation and 
yields a temperature dependence. As illustrated by the 
three curves in Fig. 6, if a ferroelectric cluster is created 
at a temperature close to T c , relaxation of this cluster is 
slow because of a rather strong dipole fluctuation, which 
holds the cluster in the metastable ferroelectric state 
from going back to the paraelectric one. Away from 
T c , t r decreases considerably for the dipole fluctuation 
is highly suppressed. This behavior is nothing but the 
critical slowing down of photorelaxation. 

In Fig. 6, it can also be seen that with the increase of 
S, T c moves to the high temperature side so as to over- 
come a higher potential barrier between the ferroelectric 
and paraelectric phases. Furthermore, the evolution of t r 
becomes gentle as well, implying a gradual weakening of 
dipole fluctuation at high temperature region. 

In Fig. 7, we show the temperature dependence of t r 
for different £ when T > T c , where S is fixed at 4.439. 
The values of parameters for this calculation are given in 
Table II. When £ changes from 3.261 to 2.861, as shown 
in Fig. 7, T c gradually increases. This is because with the 
decrease of £, the ferroelectric state at x c (refer to Fig. 3) 



becomes more stable and can survive even larger thermal 
fluctuation. In a manner analogous to Fig. 6, the evolu- 
tion of t r also displays a sharp decline at low temperature, 
and becomes more and more smooth as temperature in- 
creases. 

In Fig. 8, we plot the temperature dependence of t r 
for different barrier heights, i.e., S + £ varies from 7.0 to 
8.0, while the ratio <5/£ is fixed at 1.5. Parameters for 



a. 




T(K) 



FIG. 8: Temperature dependence of relaxation time t r for 
various barrier height 8 + £ when T > T c , where 8/^=1.5 is 
assumed. Error bars show the 95% confidence interval. 
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TABLE III: Parameters adopted for calculation of Fig. 8. 



C4 


C6 


5 




5 + t 


T c (K) 


2.1557xl0~ 2 


3.7309 xlO" 4 


4.200 


2.800 


7.000 


372 


2.0122xl(T 2 


3.2505xl0~ 4 


4.500 


3.000 


7.500 


400 


1.8860 xl(r 2 


2.8557xl0" 4 


4.800 


3.200 


8.000 


436 



this calculation are provided in Table III. As already dis- 
cussed with Figs. 6 and 7, larger S tends to increases T c , 
but higher £ applies an opposite effect on T c . Combining 
these two effects, in Fig. 8, one finds that T c increases if 
both 6 and £ are enhanced, indicating that in this case, 
the change of S plays a more significant role. Meanwhile, 
in contrast to Figs. 6 and 7, all the three curves in Fig. 
8 present smooth crossovers on decreasing temperature 
towards T c , signifying that dipole fluctuation can be pro- 
moted by lowering £ even the temperature is decreased. 

In Namikawa's experiment, the wavelength of soft x- 
ray laser is 160 A, hence the photo-created cluster is 
of nano scale. However, it should be noted that relax- 
ation of nano-sized cluster is beyond our present quantum 
Monte Carlo simulation because of the size limitation of 
our model. This is the primary reason why the exper- 
imentally measured relaxation time can reach about 30 
picoseconds, being several times longer than our calcu- 
lated results. In spite of the difference, our calculation 
has well clarified the critical dynamics of BaTiOa and the 
origin of speckle variation. 



IV. SUMMARY 

We carry out a theoretical investigation to clarify the 
dynamic property of photo-created ferroelectric cluster 
observed in the paraelectric BaTiOs as a real time cor- 
relation of speckle pattern between two soft x-ray laser 
pulses. The density matrix is calculated by a pertur- 
bative expansion up to the fourth order terms, so as to 
characterize the time dependence of scattering probabil- 
ity. The cluster-associated phonon softening as well as 
central peak effects are well reproduced in the phonon 
spectral function via a quantum Monte Carlo simulation. 
We show that the time dependence of speckle pattern is 
determined by the relaxation dynamics of photo-created 
ferroelectric cluster, which is manifested as a central peak 
in the phonon spectral function. The photorelaxation of 
ferroelectric cluster is featured by a critical slowing down 
on decreasing the temperature. Near the T c , cluster exci- 
tation is stablized by the strong dipole fluctuation, corre- 
spondingly the relaxation becomes slow. While, at higher 
temperature, dipole fluctuation is suppressed, ending up 
with a quicker relaxation of cluster. Our simulation also 
illustrates that the critical slowing down and dipole fluc- 
tuation are subject to the chemical environment of crys- 
tal. 
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